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Abstract — An interactive function computation problem in a 
collocated network is studied in a distributed block source coding 
framework. With the goal of computing a desired function at the 
sink, the source nodes exchange messages through a sequence 
of error-free broadcasts. The infinite-message minimum sum- 
rate is viewed as a functional of the joint source pmf and is 
characterized as the least element in a partially ordered family 
of functional having certain convex-geometric properties. This 
characterization leads to a family of lower bounds for the infinite- 
message minimum sum-rate and a simple optimality test for any 
achievable infinite-message sum-rate. An iterative algorithm for 
evaluating the infinite-message minimum sum-rate functional is 
proposed and is demonstrated through an example of computing 
the minimum function of three Bernoulli sources. 

I. Introduction 

In this paper, we study, using a distributed block source cod- 
ing framework, an interactive function computation problem in 
a collocated network where nodes take turns to broadcast mes- 
sages over multiple rounds. Consider a network consisting of 
m source nodes and a sink node. Each source node observes a 
discrete memoryless stationary source. The sources at different 
nodes are independent. The sink does not observe any source 
and needs to compute a samplewise function of all the sources. 
To achieve this objective, the nodes take turns to broadcast t 
messages in total. Nodes are collocated, meaning that every 
message is recovered at every node without error. After all 
the message broadcasts, the sink computes the samplewise 
function. The communication is said to be interactive if t > m. 

For all finite f, a single-letter characterization of the set of all 
feasible coding rates (the rate region) and the minimum sum- 
rate was provided in |1| using traditional information-theoretic 
techniques. This, however, does not lead to a satisfactory 
characterization of the infinite-message limit of the minimum 
sum-rate as the number of messages t tends to infinity. The 
objective of this paper is to provide a "limit-free" characteri- 
zation of the infinite-message minimum sum-rate, i.e., it does 
not involve taking a limit as t — > oo, and also an iterative 
algorithm to evaluate it. This result is similar to that provided 
in 0, where a two-terminal interactive function computation 
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problem was studied. The infinite-message minimum sum-rate 
is the fundamental limit of cooperative function computation, 
where potentially an infinite number of infinitesimal-rate mes- 
sages can be used. While the asymptotics of blocklength, rate, 
quantizer step-size, and network size have been explored in 
the distributed source coding literature, asymptotics involving 
an infinite number of messages has not, to the best of our 
knowledge, been studied and is not well understood. 

In this paper, we view the infinite-message minimum sum- 
rate as a functional of the joint source pmf. The main result 
is the characterization this functional as the least element 
in a partially ordered family of functionals having certain 
convex-geometric properties. This characterization does not 
involve taking a limit as the number of messages goes to 
infinity. The proof of this main result suggests an iterative 
algorithm for evaluating the infinite-message minimum sum- 
rate functional. We demonstrate this algorithm through an 
example of computing the minimum function of three sources. 

Related interactive computation problems in various net- 
works have been studied in |3]]-||6] using the framework of 
communication complexity Q, (8), where computation is 
required to be error-free. A function computation problem in 
a collocated network is studied in (9l within a distributed 
block source coding framework, under the assumption that 
conditioned on the desired function, the observations of source 
nodes are independent. Multiround (interactive) function com- 
putation in a two-terminal network is studied in (2), IflOll . 
IfTTII within a distributed block source coding framework. 
The impact of transmission noise on function computation is 
considered in lfl2l - lfT4l but without a block coding rate. 

The rest of this paper is organized as follows. In Sec. [II] we 
setup the problem and recap previous results. In Sec. [Til] we 
provide the main result, a "limit-free" characterization of the 
infinite-message minimum sum-rate. In Sec. [TV] we present 
an iterative algorithm for evaluating the minimum sum-rate 
functional and demonstrate it through an example. 

II. Interactive Computation in Collocated Networks 

A. Problem formulation 

Consider a network consisting of m source nodes numbered 
1, ...,m, and an un-numbered sink (node). Each source node 
observes a discrete memoryless stationary source taking values 
in a finite alphabet. The sink has no source samples. For 



each j = l,...,m, let Xj := (X/l), . . . ,Xj(n)) e (*;)" 
denote the n source samples which are available at node- 
j. In this paper, we assume sources are independent, i.e., 
for = (Xi(i),X 2 (i),...,X m (i)) are iid p x »> £ fx™ 

where Px m '•= {lYj=\ Px,] is the set of all product pmfs on 
X\ X . . . x X m . We adopt this assumption for two reasons: (1) 
to isolate the impact of the structure of the desired function 
on the efficiency of computation, (2) to obtain an exact 
characterization of the optimal efficiency. The general problem 
where the sources are dependent across nodes is open. Let 
/ : X\X. . .xX,,, — > X. be the function of interest at the sink and 
let Z(0 := f{Xy(i), . . .,X,„(i)). The tuple Z := (Z(l), . . . ,Z(n)), 
which denotes n samples of the samplewise function, needs 
be computed at the sink. 

The communication is initiated by node-A:. The nodes take 
turns to broadcast messages in f steps. In the r'-th step, node- 
j, where j - (k + i — 1 mod m), □ generates a message as 
a function of the source samples Xj and all the previous 
messages and broadcasts it. Nodes are collocated, meaning 
that every broadcasted message is recovered without error at 
every node. After f message broadcasts, the sink computes the 
samplewise function based on all the messages. If t > m, the 
communication is multi-round and will be called interactive. 

Definition 1: A f-message distributed block source code 
for function computation initiated by node-A: in a collo- 
cated network with parameters (t,n, \M\\, . . . , \Mt\) is the 
tuple (e\, . . ., e t ,g) consisting of t block encoding functions 
e\, . . . , e, and a block decoding functions g, of block-length «, 
where for every ; = 1 , . . . , f, j = (k + i - 1 mod m), 

i-l I 

ei : {Xj)" x (g) At, -> Mi, g : (g) M -» Z. 
i=] i=i 

The output of e„ denoted by M,, is called the j'-th message. 
The output of g is denoted by Z. For each ;, (l/«)log 2 |A1;| is 
called the r'-th block-coding rate (in bits per sample). 

Remark 1: (i) Each message M; could be a null message 
(|Af,| = 1). By incorporating null messages, the coding scheme 
described above subsumes all orders of messages transfers 
from m source nodes, and a f-round coding scheme subsumes 
a f'-round coding scheme if t' < t. (ii) Since the information 
available to the sink is also available to all source nodes, there 
is no advantage in terms of sum-rate to allow the sink to send 
any message, (iii) Although the problem studied in 1 1 ] is a 
special case with k — 1 and t = mr, where r e Z + is the 
number of rounds, the characterizations for the rate region and 
the minimum sum-rate in JT] naturally extend to the general 
problem described above. 

Definition 2: A rate tuple R = {R\,...,R t ) is admissi- 
ble for f-message function computation initiated by node-£ 
if, Ve > 0, 3 «(e, f) such that V« > n(e, t), there exists 
a f-message distributed block source code with parameters 
(t,n, |Ati|, ...,\M t \) satisfying 

V/ = 1, . . . , f, - log, \Mi\ < Ri + e, P(Z + Z)<€. 
n 

- j = (k mod ni) means that j s {!,,., ,m} and m divides (j ~ k). 



The set of all admissible rate tuples, denoted by 'R.,, is called 
the operational rate region for f-message function computation 
initiated by node-k. The minimum sum-rate R k sum , is given by 

mm R£Rf (Zi=i ^')- Th e f° cus °f this paper is on the minimum 
sum-rate rather than the rate region. 

Remark 2: (i) We allow the number of messages t to be 
equal to and abbreviate R k aim() to R SU m,o because there is 
no message transfer and the initial-node is irrelevant, (ii) For 
t < m, function computation may be infeasible, i.e., may 
be empty. If so, we define R k sumt := +°°. For special px»< and 
/, however, computation may be feasible even with t < m; 
in that case, R k sumt would be finite, (iii) For all t e Z + , 
^s»m r — Rsum r+T - holds, because the last r messages could 
be null. Hence the limit lim,^ M Br , =: B k sumoo exists and is 
finite, (iv) For all r E Z + , R k mm t > fl*7,v ™ d holds, because 
the first r messages could be null. It follows that R k nm x is 
independent of k and we abbreviate it to /? s „„,,co. For all finite f, 
however, we keep the superscript in R k um t because this notation 
is convenient in the proof of Theorem [T] 

For all finite f, a single-letter characterization of "R^ and 
R k umt was provided in Theorem 1 and Corollary 1 of |[T|| . 
This, however, does not directly lead to a satisfactory char- 
acterization of the infinite-message limit R s „ m ,ca, which is a 
new dimension for asymptotic-analysis involving potentially 
an infinite number of infinitesimal-rate messages. The main 
contribution of this paper is a novel convex-geometric charac- 
terization Of R S um,ca- 

B. Characterization of R k uml for finite t 

Fact 1: (Characterization of R k sumt Hj] Corollary 1]) 

R k sumt = min I(X'";U'), (2.1) 

where f^ipx-) is the set of all pu'\x"> such that (i) 
H(f(X m )\U<) = 0, (ii) V/ E [l,...,t},j = (k + i-l 
mod m), Ui - (f/ i_1 ,Z^ - (X^X'J^), and (iii) the cardinalities 
of the alphabets of the auxiliary random variables U' are 
upper-bounded by functions of |<Yi|, . . . , \X m \ and f. 

The Markov chain conditions in Fact Q] are equivalent to the 
following factorization of pu'\x m '- 

Pu'\x<» = Pu,\x t ■ Pu 2 \x tM mod ,„,£/, ■ Pu 3 \x it+2 mol]m ,u 2 ■■■■ (2-2) 

The cardinality bounds in Fact Q] which can be derived 
using the Caratheodory theorem are omitted here for clarity. 
Although the exact expressions of the cardinality bounds are 
unimportant for our discussion, a key property that needs to be 
highlighted is that the bound on the alphabet of U, increases 
exponentially with respect to (w.r.t.) f. Therefore the dimension 
of the optimization problem in 12. 11 explodes as f increases. 

Using Fact [1] we could compute R k um t for a large f to ap- 
proximate Rswn.oz- This is impractical because (i) the dimension 
of the optimization problem is large, (ii) the characterization 
of R k sum t does not inform us how close R k llm , is to R sum ,co. 
Alternatively, we could compute R k um t for increasing values 
of f until \R k , - R k ,1 falls below a threshold. However, 

1 sum,t— 1 sum,t> ' 



the dimensionality of the optimization problem grows expo- 
nentially with increasing values of t and there is no obvious 
way to reuse the computations done for evaluating R k umt _ l 
when evaluating R k sum ,. Finally, if we need to evaluate R SU m,°o 
for a different px>", we need to repeat the entire process. 

In Sec. [Ill] we take a new fundamentally different approach. 
We first view R SU m,oa as a functional of px m for a fixed /. 
Then we develop a convex-geometric blocklength-free char- 
acterization of the entire functional R SU m,°o(Px m ) which does 
not involve taking a limit as t — > oo. This leads to a simple 
test for checking if a given achievable sum-rate functional 
of px<» coincides with R S i, m ,co(px : ")- It also provides a whole 
new family of lower bounds for R slm ,<x>- In Sec. |IV] we use 
the new characterization to develop an iterative algorithm 
for computing the functional R sum ,oo(px'») and R k unh ,{px-) (for 
any finite t) in which, crudely speaking, the complexity of 
computation in each iteration does not grow with iteration 
number, and results from the previous iteration are reused 
in the following one. We demonstrate the iterative algorithm 
through an example. 

III. Characterization of R s „ m ,<x,(px m ) 
A. The rate reduction functional p k (px m ) 

If the goal is to losslessly reproduce the sources, the min- 
imum sum-rate is equal to H(X m ) = Y!k=\ H(Xk) because the 
sources are independent. The minimum sum-rate for function 
computation cannot be larger than that for lossless source 
reproduction. The reduction in the minimum sum-rate for 
function computation in comparison to source reproduction 
is given by 



p* :=H(X m )-R k 



f - max H(X m \U'). (3.3) 



A quantity which plays a key role in the characterization of 
Rsum,°o is po - the "rate reduction" for zero messages (there 
are no auxiliary random variables in this case). Let 

Pf.= {p x »ePx»:H(f(X m )) = 0}. 

Error-free computations can be performed without any mes- 
sage transfers if, and only if, px<<< e Pf- Thus, 



R. 



Po 



i,0 - 



0, if p x » e P f , 
+oo, otherwise, 



H(X m ), 



if Px» e P f , 
otherwise. 



(3.4) 



Remark 3: If /(x m ) is not constant, for all px»< e Pf, we 
have supp(/?x»') #^iX...x X m . Such px<« can only lie on the 
boundary of Px<*- 

Evaluating R k umJ is equivalent to evaluating the rate re- 
duction p k t . It turns out, however, that p^ :- lim^oopf = 
H(X"') - R sllm ,co is easier to characterize than R sum ,oo (see 
Remark ??). The rate reduction functional is the key to the 
characterization. 



B. Main result 

Generally speaking, p k , po, and p^ depend on px<<< and 
/. We will fix / and view p k (p x -), Po(px»), and pUpx™) 
as functionals of px<<< to emphasize the dependence of px*. 
Instead of evaluating poo{px m ) for one particular px<» as it is 
done in the numerical evaluation of single-terminal and Wyner- 
Ziv rate-distortion functions, our approach is to characterize 
and evaluate the functional poa(px m ) for the entire set of 
product distributions Px« rather than for one particular px<«- 
To describe the characterization of the functional pco(px m ), it 
is convenient to define the following family of functionals. 

Definition 3: (Marginal-distributions-concave, po-major- 
izing family of functionals T) The set of marginal- 
distributions-concave, po-majorizing family of functionals T 
is the set of all the functionals p : Px m — > R satisfying the 
following conditions: 

1) po-majorization: Vp x » e Px-, p(px»0 > Po(px-)- 

2) Concavity w.r.t. marginal distributions: For all k e 
{1, . . . ,m}, with p X j held fixed for all j + k,p (l~["Li Px) 
is a concave function of px k - 

Remark 4: Since po(px») = -°° for all px<« t Pf, condition 
1) of Definition [3] is trivially satisfied for all px<» e Px m \ 
Pf (we use the convention that Va e R, a > -oo). Thus the 
statement that p majorizes po on the set Px<« is equivalent to 
the statement that p majorizes H(X' n ) on the set Pf. 

Remark 5: Condition 2) does not imply that p(px m ) is 
concave w.r.t. the joint pmf px<»- In fact, Px<<< is not convex. 

We now state and prove the main result of this paper. 

Theorem 1: (i) p M e T ' . (ii) For all p e f, and all px<« € 
Pxm, we have pco(px»0 < p{px"). 

The set T is partially ordered w.r.t. majorization. TheoremQ] 
says that p M is the least element of f ' . Note that there is 
no parameter t which needs to be sent to infinity in this 
characterization of p x . 

To prove Theorem[TJwe will establish a connection between 
the f-message interactive coding problem and a (f- l)-message 
subproblem. Intuitively, to construct a f-message interactive 
code with initial-node k and px>» = Y\T=i Px t , we need to begin 
by choosing the first message which corresponds to choosing 
the auxiliary random variable U\. Then for each realization 
U\ = «i, constructing the remaining part of the code becomes 
a (t - l)-message subproblem with initial-node k + := (k + 1 
mod m) with the same desired function, but with a different 
joint source pmf px»\u, = n!=i Px > where for all i + k, p' x = 
pxi and p' x = px k \Ui ■ We can repeat this procedure recursively 
to construct a (t— l)-message interactive code. After t steps of 
recursion, we will be left with the trivial 0-message problem. 

Proof: (i) We need to verify that p ro satisfies the two 
conditions in Definition [3] 

1) Since Vpx»' 6 Px™, R sum. Co {px») < R.mmfi{px"<), we have 
Pcc(px-) > po(px»0- 

2) For any k e {l,...,m}, consider two arbitrary distribu- 
tions px h o and px k ,\, and arbitrary distributions px s for all 
j + k. For mi = 0,1, let px» Ml := p Xt , Ul ■ UJ=i,j^Px r For 



A e (0, 1), let px«<,A '■- Apx»>,\ + (1 — X)p x ™fl- We will show that 
Poo(px»;a) > Apoo(px»',\) + (1 - X) Poo(px*,o)- Let U\ ~ Ber(A) 
and (X'", t/*) ~ px»\ Ui Pu\(u\), which imply p x * = px«\x e fx™ 
and px»<\u"(-\u\) = px™, Ui e r Px<»- For all t e Z + we have, 



P?(Px»vO 



max Z/(X" ! |t/') 



= max <; 

Pu,\x k 



Co) 



max H(X"'\U') 



> max Z/(X'"|t4 

/>^ix f />u<ix»<u; ep fO>A™..i) 

= A- max H(X m \U' 2 ,U; = 1) 

/'t/;ix 4 P^ix'» I /; e ' p f(/'x»M) 
+ (1 - /}) ■ max H(X m \U' 2 , U\ = 0) 

Pj4PP"(7«(-|-.0): 

/'f/;ix t /'^ix™(/; e!P f(w".o) 

( => Ap^Cp^O + Cl-^p^C^o). (3.5) 

In step (a) we replaced pui\x k with the particular pu\\x k defined 
above. Step (b) follows from the "law of total conditional 
entropy" with the additional observations that conditioned 
on U* = «i, px»|£/!f(-|"i) = px", m and i/(X m |t/^, f/* = u x ) 
only depends on Pt/^|X'»(/;( p h Step (c) is due to the ob- 
servation that for a fixed pu*\x r conditioned on U* { = u\, 

Pu\\x k Pu 2 \x>»u\ 6 ^(Px-, Ul ) iff P(/;n»(/; e Pf^ipx-,,,,)- Now 
send r to infinity in both the left and right sides of (13.51 . 
Since lim^oopf = lim, 

— >co — Poo-, 

we have Poo(pxv) > 
^Pco(px'»,i) + (l-A)pco(/?x»',o)- Therefore, poo satisfies condition 
2) in Definition [3] Thus, p m e !F. 

(ii) It is sufficient to show that: Vp e T, Vpx»< e fx™, 
St e Z + UiO}, V& € {1, . . .,m}, p*(px») < p(px»). We prove this 
by induction on t. For f = 0, the result is true by condition 1) in 
Definition [3] Assume that for an arbitrary t e Z + , pf^px") - 
p(px»0 holds. We will show that p k (px m ) < p(px m ) holds. 

p*(px») = max #(X m |[/') 



max<! max H(X"'\U') 

PUitx t Pu' 2 \x'«u^ 

Pu 1 \x t Pu' 1 \x>»v i e '? J :(Px>») 



id) 



(e) 



max { 



^ P£/i("l) 



max //(X m |t/^, t/i = Mj) 

Pu'^x'-ut 



max i y, Ptfi("l)PM(PX»|tfi01"l)) 
^i,e supptpuj) J 

< max] V PuS.m) p{px»<\uM u ^)\ 

Pu,\x k ' 



(3.6) 



max I y pt/, («i) p(Pxm ( - l M i) Px'-'XE. , ) \ 
P u \i x k , > 

^Hiesupp(poj) J 



(A) 

< max •! p 

Pv,\x k 
= P(PX"0- 



^ puMdPxtW^ui) Px<-<x»;, 

(,«l€ SlippOyj) 



The reasoning for steps (d) and (e) are similar to those for 
steps (b) and (c) respectively in the proof of part (i) (see 
equation array d3.51 l). In step (e) we need to use the fact that 
Px»|£/i('|wi) £ Px'», which is due to i2.2\ and the assumption 
that px>» € Px m - Step (f) is due to the induction hypothesis 
p k t _ y (px"<) ^ p(px m ) for all k. Step (g) is due to the Markov 
chain U\ - X k - {X k - l X" k ' +y ) and because X k and (X*" 1 ^) 
are independent. Step (h) is Jensen's inequality applied to 
P(Px k • Px'-'x™,) which is concave w.r.t. px k - ■ 
Since every pe f gives an upper bound for p M , (H(X m )-p) 
gives a lower bound for R sum ,<=o- This fact provides a method 
for testing if an achievable sum-rate functional is optimal. If 
R*(px») is an achievable sum-rate functional then Vpx» e 
R\px»<) ^ R S um,co(Px>»)- If it can be verified that p* := (H(X m )- 
R*) e f, then by Theorem [TJ R* = R sllmtC o- 

IV. Iterative algorithm 

Although Theorem [TJ provides a characterization of p m and 
Rsum,co that is not obtained by taking a limit, it does not 
directly provide an algorithm to evaluate R sum ,<x- To efficiently 
represent and search for the least element of T is nontrivial 
because each element is a functional; not a scalar. The proof 
of Theorem [TJ however, inspires an iterative algorithm for 
evaluating R k sum t and J?,„ m>0O . 

Equation ( 13.61 ) states that p*(px m ) is the maximum value 
of p e R such that (px'»,p) is a finite convex combination 
of {(px^Uii-lmXp^ipxm^i-M)^ sup V (p Vl ), where p x .»(-) and 
Px"<\Ui('\ u i) have the same marginal distributions px for all j + 
k and differ only on px k ■ Now we fix the marginal distributions 
px for all j + k, and consider the hypograph of p*_j w.r.t. px^ 
hyP-pf-i := i(Px k ,p) ■ P < p£i(Tl? =i Px,)}. Due to the 



convex hull of hyp 



PxA-l 



is hyp p t . This relation enables 



us to evaluate p, from p l { : fixing px ; for all j + k, p. 



is 



the least concave functional w.r.t. p Xk that majorizes pj _ r In 
the convex optimization literature, (— p*) is called the double 
Legendre-Fenchel transform or convex biconjugate of (— p*_j) 
|fT31 . We have the following iterative algorithm. 
Algorithm to evaluate R k sum>t 

. Initialization: For all k = l,...,m, define p k {px m ) = 

Po(Px»0 by equation (E§ for all p x »- in V x - = {n"=i Px,}- 
• Loop: For t = 1 through f do the following. 

For every k = 1, . . . ,m do the following. 

For every set of marginal distributions {px,}™-! jrf do the 

following. 

- Construct hyp^p*^. 

- Let p\ be the upper boundary of the convex hull of 

^Px/r-V 



. Output: R k swn>t (px») = H{X m ) - p k t { Px ,„). 

To make numerical computation feasible, Vx™ has to be 
discretized. Once discretized, however, in each iteration, the 
amount of computation is the same and is fixed by the dis- 
cretization step-size. Also note that results from each iteration 
are reused in the following one. Therefore, for large t, the 
complexity to compute R sum , grows linearly w.r.t. t. 

R sum ,oo can also be evaluated to any precision, in principle, 
by running this iterative algorithm for t = 1,2,..., until 
some stopping criterion is met, e.g., the maximum difference 
between p k { and p k on Px m falls below some threshold. De- 
veloping stopping criteria with precision guarantees requires 
some knowledge of the rate of convergence which is not 
established in this paper and will be explored in future work. 
When the objective is to evaluate R SU m,oo(px»>) for all pmfs 
in fx m , this iterative algorithm is much more efficient than 
using d2.lt to solve for R k sumt for each p x <« for t = 1,2,..., 
an approach which follows the definition of R su „ ha3 literally as 
the limit of R k umt as t — > oo. Our iterative algorithm is based 
on Theorem [T] which is a characterization of R SLtm ,oo without 
taking a limit involving t. 

Example: (MIN function) Take m = 3 nodes. Xj ~ Ber(p,). 
/(x 3 ) = min, = i 2,3 x,-. The joint pmf px> is parameterized by 
p = (pi,p2,Pi) £ [0, l] 3 . It is easy to see that 



Fig. [2] shows the plots of the rate reduction function with 



Po(p) 



ZLi h 2 (pd, 



ifpe P f , 
otherwise, 



where P f = {p : p x = p 2 = py = 1 or pip 2 p3 = 0}. 

Now let us fix p Xl and p Xl and apply the convex biconjugate 
operation on po w.r.t. p X} to obtain p\. Specifically, for every 
fixed (p\,p2), we focus on po on the line segment {pi}x{p2}x 
[0, 1] and convexify hyp ft p = {(p3,p) ■ p < Po(pi,P2,pi)} to 
obtain hyp ft pj. Then we repeat this procedure but applying 
the convexification operation w.r.t. px 2 , px,, etc to obtain p\, 
pi, etc. In numerical computation, p takes values on a discrete 
grid where p\,pi,pz are multiples of a finite step size A. The 
convexification operation involves finding a convex hull of a 
finite number of points in a plane. 

As we decrease A and increase t, p, approximates p M . 
Fig. Q] shows p t (l/2, 1/2, 1/2) for different t and A. For each 
A, p t converges as t increases. For a small enough A (fine 
enough discretization), the limit represents the actual value 
of poo(l/2, 1/2, 1/2). Notice that for a small enough A, p t 
keeps increasing as t grows, which means there is always an 
improvement for using more messages. 
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Fig. 2. Rate reduction function for t = 40 
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Fig. 1. p,(l/2, 1/2, 1/2) for different step sizes A 



